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In many cases of interest such as a multi- 
layer mesoscopic structure (e.g. a superlattice) 
or the grain boundaries network in high-T c 's, 
one deals with the situation where electrons tra- 
verses many partially transparent interfaces with- 
out loosing coherence. Then, the consequent re- 
flection/transmission events require a simultaneous 
consideration. Theoretically, even an isolated inter- 
face poses certain difficulties: Since abrupt changes 
violate the quasiclassical condition, the theory of 
superconductivity in terms of the quasiclassical ma- 
trix Green's function g R |IJ is invalid at interfaces. 
The interface is included via the boundary condi- 
tion derived by Zaitsev H — a cubic matrix relation. 
In case of many interfaces, one comes to a system of 
nonlinear matrix equations solution of which is far 
not simple if possible. Moreover, some authors ar- 
gue P] that the normalization condition is violated 
in the multiple interface case so that the quasiclas- 
sical scheme fails. The purpose of this paper is to 
suggest a new method which allows one to study 
systems with multiple partially transparent inter- 
faces SIS' or SIN in the framework of the quasi- 
classical theory of superconductivity. Motivated by 
the paramagnetic effect problem (see references in 
H ) , the response to the vector potential (the super- 
fluid density) of a I SIS' I sandwich is calculated. 

A convenient starting point is the formulation 
|0| of the quasiclassical technique in terms of the 2- 
point trajectory Green's function g R (xi,x 2 ) where 
xi^ are the coordinates along the trajectory char- 
acterised its direction n and position R _L n. The 
equation of motion reads 

{iv—+ H R (x 1 ))g R (x 1 ,x 2 ) =iv5{x! - x 2 ) (1) 
where H R = 



A 



— Y, R and the self-energy is 

denoted ~Z R . The Green's function g R Fn of Ref. @] 
is found by g R (x ± 0, x) = (g R Fn ± l)/2. 

The key point is the understanding that in the 
presence of partially reflective interfaces a typical 
trajectory remains extended and simply connected. 
In other words, the trajectories typically do not 
have loops and are not split to interfering paths 
so that there is only one path connecting any two 
points on the trajectory. Then, the coordinate along 
the trajectory is well defined and the method of 
2-point Green's function is applicable. Due to re- 
flection/transmission at the interfaces, the typical 
path becomes tree-like with "knots" , i. e. the points 
where ballistic pieces of the trajectory cross. In 
between knots, the quasiclassical equation Eq.(Q) 



is valid, whereas the knots are included via the 
boundary matching condition. 

The latter is formulated for the "wave functions" 
= 0' = ( u > ~ u ) factorizing g R , 

g R (x + O,x)=0+(x)^(x), 

found from the linear equations, 



(iv 



^+H R {x))<f> = Q ,<t>±(x^±oc)=0, (2) 

and normalized, = 1. 

A general interface is partially diffusive, mixing 
states within a continuum of momenta. As a model, 
one allows mixing a finite number of states, in- 
coming and N outgoing, at a knot. The knot value 
of the wave functions <f>i (4>k' ) in the in-coming (out- 
going) channels, i, k' = 1, . . . , N are related to each 
other by the unitary scattering matrix [^) 
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(3) 



The energy independent matrix S, the normal state 
property, is considered as a given input. 

In this scheme, one finds the "wave function" (j> 
by solving Eq.(||) on the pieces of the trajectory 
in between the knots and uses Eq.(||) to tailor the 
pieces. The scheme allows one to consider the usual 
and Andreev reflections on equal footing, simul- 
taneously with (in)elastic scattering processes in- 
cluded via the self-energy T, R . 

Here we restrict ourselves to the simplest case of 
N = 2 (corresponding 
e.g. to a partially trans- 
parent specular interface) 
when the S-matrix takes 
the form ( JJ, r *. V R = 

M 2 , T = \t\ 2 , R + T = 1. 
Analysis Q shows that Eq.(|J) leads to the follow- 
ing relation among the knot values of the 1-point 
Green's functions g R and g R on the trajectories 1 
and 1': 

g^MgfM- 1 (4) 
where the transfer matrix M. reads 
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here g R %2 i the "across the interface" propagator, is 



(6) 



where [...]( + )_ denotes (anti)commutator. 



Eqs.(|||3|) being equivalent to Zaitsev's boundary 
conditions, seem to be more convenient for numer- 
ical calculations. Besides, they can be easily gener- 
alised || to the case of a diffusive interface mod- 
elled by a knot N > 2. 

To illustrate the method, we calculate the den- 
sity of states and the superfluid density for a sys- 
tem with 3 planes of reflection: a sandwich made of 
two superconducting layers L and R separated by 
a partially transparent interface; the layers occupy 
the regions — d^ < z < and < z < d R and the 
order parameter is Al and Ar, respectively. 

Consider the periodic zig-zag trajectory in the R- 
layer, formed by reflections from the boundary at 
z = dfi and the interface at z = 0; it is specified by 
the angle and the period ae = 2d R / cos8. Knots 
are the points where the zig-zag touches the inter- 
face. Because of the periodicity, each of the regu- 
lar solutions <f)±(x) is proportional to one of the 
eigenvectors of the evolution matrix U x defined via 
<p(x + ae) = U x 4>{x). The evolution matrix reads, 

U x = U { ° ) (x,0)MU ( - o) (0,x-a e ) , (7) 

where M. is the transfer matrix Eq.(||), and 
is the bulk evolution operator corresponding to 
Eq.(^): U(°\x+y,y) — exp[ixH R /v]. One can show 
H that the 1-point Green's function for the direc- 
tion n z = ± cos 9 can be calculated as 



g R (x) = U' x 
where UL stands for U, 




(8) 



Given the interface value of the Green's function 
in the L-layer, one finds the transfer matrix from 
Eqs.(||j6|), and Eqs(Q,||) give the R-layer Green's 
function. Equations analogous to can be writ- 
ten for the L-layer. Numerical iteration of the sys- 
tem of the equations, allows one to evaluate the 
physical properties of the system. 

Motivated by the recent idea || about the 
paramagnetic instability at normal-metal - super- 
conductor interfaces in the situation where the 
proximity induced order parameter in the normal 
metal is negative (repulsion), we consider the case 
when Al — — Ar (questions related to the self- 
consistency equation are left aside here). The angu- 
lar resolved density of states (DOS), i.e. Re (g R )u 
for the almost transparent interface R = 0.1 is pre- 
sented in Fig. 1. Note the zero value of DOS at 
small energies: DOS is extremely sensitive to reflec- 
tion which leads to the splitting of the zero energy 
Andreev levels. Accordingly, the superfluid density 
p s (T) (j = -const- p s {T)A), 

Ps = l-J-^ (-^) (MS R )n) 
((...) denoting the angular and space averaging) 
is also very sensitive to reflection. As seen from 
Fig. 2, the low temperature paramagnetic effect (i.e. 
p s < 0) being strong for ideal interface, R = 0, dis- 
appears when the probability of reflection is as low 



as 0.04. These findings confirm the theoretical pos- 
sibility of the paramagnetic effect, and at the same 
time indicate that the effect is very sensitive to the 
inteface reflection inevitable in the experiment. 

In conclusion, a new method which allows one 
to describe a system with multiple interfaces in the 
framework of the quasiclassical theory is presented. 




Figure 1: Angle resolved DOS vs (energy/A); 
A L = -A R = A, d L /cosd = d R /cos6 = v/A 
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Figure 2: Space averaged superfluid density vs tem- 
perature for the reflection R = 0, 0.01, 0.04 
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